Directional (version 3.1)

Random values generation from the ESAG distribution: Random values generation from the ESAG distribution

Description

Random values generation from the ESAG distribution.

Usage

ESAGsim(n, param)

Arguments

n

A number; how many vectors you want to generate.

param

The parameters of the ESAG distribution, the first three are the mean vector in $R^3$ and the next two are the two gammas.

Value

An nx3 matrix with the simualted unit vectors.

Details

A random sample from the ESAG distribution is generated. In case the gammas are zeo the sample is drawn from Independent Angular Gaussian (IAG) or projected normal.

References

Mardia, K. V. and Jupp, P. E. (2000). Directional statistics. Chicester: John Wiley & Sons.

Paine P.J., Preston S.P., Tsagris M. and Wood A.T.A. (2017). An Elliptically Symmetric Angular Gaussian Distribution. Statistics and Computing (Accepted).

See Also

ESAGmle, ESAGdensity, spml.mle, acg, circ.summary

Examples

# NOT RUN {
m <- colMeans( as.matrix( iris[,1:3] ) )
y <- ESAGsim(1000, c(m, 1, 0.5) )
ESAGmle(y)
# }